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CT APPARATUS, CT IMAGING METHOD AND METHOD OF 
PROVIDING SERVICE USING THE SAME 
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BACKGROUND OF THE INVENTION 
5 The present invention relates to an 

industrial computed tomographic (CT) apparatus, a CT 
imaging method and a method of providing service using 
the same. 

In industrial X-ray CT apparatuses it is 
10 common practice to place an object of concern on a 

turntable that rotates in a horizontal plane and take a 
horizontally sliced tomographic image. When an object 
is put on the turntable of the X-ray CT apparatus, the 
user selects reference planes (horizontal plane, 
15 vertical plane, etc.) and aligns a coordinate system of 
the object with that of the turntable by using a square 
and level- In inclining the object, the user makes an 
adjustment by inserting a wedge-like jig under the 
object . 

20 An industrial X-ray CT apparatus with a 

setting device that allows an object to be set in 
different positions and directions is disclosed in JP- 
A-'2000-298106. 

In a medical X-ray CT apparatus, a method is 

25 described in JP-A-2000-107161 which involves imaging a 
patient body with markings on his or her body surface 
and side surfaces and determining a position of a 
legion as a relative position from the body surface 
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markings. 

With JP-A-2000-298106, although the 
installation of an object to be imaged is made easier 
for the user, the work for determining reference planes 
5 of the object is difficult to achieve. 

The method of JP-A-2000-107161 can be applied 
to medical apparatus because the position of a human 
body attached with markings on its surface is fixed. 
However, in industrial apparatus the reference position 
10 (reference plane) changes from one object to another, 
so this method cannot be applied . 

SUMMARY OF THE INVENTION 

An object of the present invention is to 
provide a CT apparatus, a CT imaging method and a 

15 method of providing service using these, in which, even 
when objects with different reference planes are set on 
a turntable for imaging, the user can easily set the 
objects without requiring a process of determining 
reference planes of the objects. 

20 The present invention determines a geometric 

feature (information representing a geometry) of a 
contact surface between a marker and an object to be 
imaged by using a tomographic image from a CT apparatus 
of the object attached with the marker on a reference 

25 plane; determines from the geometric feature a 

positional relation between a reference coordinate 
system of the object and a coordinate system of the CT 
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apparatus; and obtains from the positional relation a 
tomographic image (bit-map data) of a desired portion 
of the object- 

Other objects, features and advantages of the 
5 invention will become apparent from the following 

description of the embodiments of the invention taken 
in conjunction with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a diagram showing a first 

10 embodiment concerning a CT imaging method of the 
present invention . 

Figs. 2A-2C are schematic diagrams showing a 
method of attaching markers to an object to be imaged. 
Fig. 2A representing an example of markers attached to 

15 three orthogonal planes. Fig. 2B representing an 
example of markers attached to a plane and a 
cylindrical surface, crossing each other at right 
angles, and Fig. 2C representing an example of markers 
attached to orthogonal two planes and one vertex. 

20 Fig. 3 is a schematic diagram showing an 

example of an X-ray CT apparatus of this invention. 

Fig. 4 is a diagram showing an example screen 
displayed on a marker arrangement pattern presenting 
means A6. 

25 Fig. 5 is an explanatory diagram showing a 

marker identifying process P3 of Fig. 1- 

Figs- 6A-6D are explanatory diagrams showing 
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an example of how geometric features are extracted by 
the surface domain expansion method. Fig. 6A 
representing a group addition process. Fig. 6B a group 
boundary determination process. Fig. 6C a group 
5 expansion process, and Fig. 6D a state in which a group 
representative surface and boundaries are recalculated 
by the group boundary determination process. 

Fig. 7 is a diagram showing a typical example 
of a marker surface feature domain list D3. 
10 Fig. 8 is an explanatory diagram showing a 

coordinate system determination process P4 of Fig. 1. 

Fig. 9 is a schematic diagram of a second 
embodiment of the invention showing a marker attached 
to an object- 

15 Fig. 10 is a diagram showing the second 

embodiment concerning a CT imaging method of the 
invention - 

Fig. 11 is an explanatory diagram of a marker 
identifying process P3a of Fig. 10. 
20 Fig. 12 is a schematic diagram showing a 

third embodiment of the invention concerning a CT 
imaging service method. 

DETAILED DESCRIPTION OF THE EMBODIMENTS 
(First Embodiment) 
25 A first embodiment of the present invention 

will be described by referring to the accompanying 
drawings. Fig. 1 illustrates a method of taking a CT 
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image according to the first embodiment. This CT 
imaging method has a marker application process PI, a 
preliminary imaging process P2, a marker identifying 
process P3, a coordinate system determination process 
5 P4, and a full imaging process P5. 

In the marker application process PI, a 
marker is attached to a reference plane of an object to 
be imaged T. For example, when the object has a plane 
as a horizontal or vertical reference, the marker is 

10 attached to that plane. If the object has a 
cylindrical plane as a horizontal or vertical 
reference, the marker is attached to that cylindrical 
plane- The marker is made of a material with an X-ray 
transmissivity (or attenuation rate) different from 

15 that of the object. For example, a material with a 
different density from that of the object has a 
different X-ray transmissivity. The marker size 
(length and width) is preferably in a range of 20 mm to 
50 mm. The marker is also preferably made of a plastic 

20 material for easy bonding. The method of bonding the 
marker will be described later with reference to Fig. 
2. 

In the preliminary imaging process P2, the 
object T is placed on a turntable of an X-ray CT 
25 apparatus and is wholly imaged at intervals of five to 
20 times of a target tomographic pitch . For example, 
if the target tomographic pitch (final desired 
tomographic pitch) is 0.4 mm, the image is performed at 
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a greater interval of 2-8 mm. 

Once the object is imaged by the X-ray CT 
apparatus, tomographic images (two-dimensional bit-map 
data) are produced for individual cross-sections. A 
5 number of tomographic images thus obtained are stacked 
together to produce a tomographic image set (three- 
dimensional bit-map data) . In the preliminary imaging 
process P2, a coarse tomographic image set Dl, a set of 
tomograms with a coarse (low) resolution, is obtained. 

10 In the preliminary imaging process P2, to 

prevent the object from being shifted in position by 
the rotation of the turntable, the object need only be 
fixed and there is no necessity for adjusting a 
placement angle of the object. The construction of the 

15 X-ray CT apparatus used in P2 will be described by 
referring to Fig. 3. 

In the marker identifying process P3, the 
coarse tomographic image set Dl obtained in the 
preliminary imaging process P2 is subjected to three- 

20 dimensional image processing to identify the marker. 
For example, if the object being imaged is made of 
aluminum, as in the case of an engine block, its 
density is about 2.7 g/cm^. When a clay is used as the 
marker, its density is about 1.5 g/cm^. In this case, 

25 the X-ray transmissivity of the marker is greater than 
that of the object and the X-ray attenuation rate of 
the marker is smaller than that of the object. 

Since a density distribution in the interior 




of the object can be known by using the X-ray CT 
apparatus, it is possible to distinguish only the clay 
(marker) portion from the aluminum and air portions for 
its extraction- In this way, the marker is identified 
5 by taking advantage of these differences between the 
marker and the object in the density. X-ray 
transmissivity and X-ray attenuation rate. 

The clay marker is plastically deformed when 
pressed against the surface of the object. When the. 

10 surface of the object is flat, the marker's contact 
surface with the object is also flat. By performing 
plane extraction processing (processing to 
automatically recognize and extract a plane of more 
than a predetermined area) on this contact surface, 

15 geometric definitions (geometric features) of the 

contact surface can be calculated with the turntable as 
a reference. The geometric definitions of all contact 
surfaces thus calculated are recorded in a marker 
surface feature domain list D3- The plane extraction 

20 processing may employ a plane identification method 
described later, such as a surface domain expansion 
method. The marker identifying process P3 will be 
detailed by referring to Fig, 5. 

In the coordinate system determination 

25 process P4, a coordinate system D4 of the object with 
the turntable taken as a reference is calculated from 
the geometric definitions of the contact surface 
between the marker and the object, which are recorded 
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in the marker surface feature domain list D3. If a 
plurality of markers are all attached to planar 
surfaces, the individual contact surfaces are compared 
with a horizontal reference plane (XY plane) and 
5 vertical reference planes (YZ plane and ZX plane) with 
the turntable as a reference, and then an affine 
transformation matrix is determined which makes an 
overall angular difference minimum. 

For example, the affine transformation matrix 

10 is determined so that a sum of the squares of angles 
between each marker and the reference planes is 
minimum. The affine transformation matrix is a matrix 
that represents a relation between two coordinate 
systems based on a linear, combination of basic vectors. 

15 Here, matrices corresponding to a rotation and a 

translation of Cartesian coordinate systems are used. 
The detail of this process will be described with 
reference to Fig. 8. 

In the real or final imaging process P5, the 

20 object T is imaged at a target tomographic pitch to 

produce high resolution tomographic images D5 and/or a 
tomographic image set D6, a final result of the 
process- The tomographic images D5 and the tomographic 
image set D6 are bit-map data. The user (operator) 

25 inputs into the X-ray CT apparatus beforehand 

coordinate values of a portion the user wants to image. 
Since the coordinate svst em D4 of the object with the 
turntable taken as a reference is already calculated in 
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the coordinate system determination process P4, the X- 
ray CT apparatus automatically corrects, before taking 
images, the coordinate values of the target portion 
entered using the coordinate system D4 . This 
5 correction is made to determine the imaging position 
and the imaging range specified by the user in advance. 

The X-ray CT apparatus also corrects the 
number of images. This correction is made, for 
example, when the imaging cross-sections are inclined 

10 to the horizontal plane as a result of performing the 
affine transformation on the imaging cross-sections 
specified by the user. The number of images to be 
taken is determined by dividing a length in the 
vertical direction of a domain from a lowest point to a 

15 highest point of the imaging- cross-sections by a 
vertical resolution characteristic of the X-ray CT 
apparatus. In this way, based on the positional 
relation between the reference coordinate system of the 
object and the coordinate system of the X-ray CT 

20 apparatus, bit-map data of a desired portion of the 
object is obtained. 

This CT imaging method has the above- 
described processes and thus can significantly automate 
the step of determining the coordinate systems D4 of 

25 objects to be imaged which have different (arbitrary) 

reference planes. That is, this method allows the user 
to omit the step of determining the reference planes of 
the objects, thus facilitating the installation of the 
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objects and the operation of the X-ray CT apparatus. 

Examples of attaching markers to objects are 
schematically shown in Fig. 2A to Fig. 2C. There are 
different patterns of attaching the markers. These 
5 marker arrangement patterns are such as will allow the 
coordinate system determination process P4 to determine 
the affine transformation matrix- 
Fig. 2A represents a case in which markers M 
are attached one to each of three orthogonal planes of 

10 an object T. In this case, an intersecting point S of 
the three planes constitutes a reference point of the 
coordinate system and intersecting lines L1-L3 between 
two planes constitute reference axes. In this example, 
all three markers M have different sizes. This allows 

15 the marker identifying process P3 to automatically 
identify individual markers M. For the purpose of 
distinguishing between the markers, the materials of 
the markers M may be dif f erentiated- 

Fig. 2B represents a case where markers M are 

20 bonded to a flat surface Fl and a cylindrical surface 

F2 of an object T, these planes intersecting each other 
at right angles. In this case, a center axis of the 
cylindrical surface F2 forms a reference axis of the 
coordinate system, and an intersecting point of the 

25 reference axis and the flat surface Fl constitutes a 
reference point of the coordinate system. With the 
surface domain expansion method described later, since 
a flat surface, cylindrical surface, spherical surface, 
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revolution surface, quadratic surface and pushed or 
translated surface can be recognized and extracted, the 
coordinate system can be defined by using these 
surfaces- The revolution surface is a locus obtained 
5 by rotating a line (straight line, curved line, and a 
combination of these) about one straight line (center 
axis) . The quadratic surface is a surface described by 
a quadratic equation of coordinate variables (X, Y, Z) . 
The translated surface is a plane of locus obtained 

10 when a line is moved parallelly. 

Fig. 2C represents a case where markers M are 
bonded to two orthogonal flat surfaces of an object T 
and to one vertex. In this case the vertex attached 
with the marker M constitutes a reference point of the 

15 coordinate system, and an intersecting line LI between 
the two planes constitutes a reference axis. If the 
markers M are bonded to a vertex and an edge line of an 
object T, the contact surface between the markers M 
and the object T is not on one plane but on three or 

20 two divided planes. By identifying and extracting all 
of these contact surfaces in the marker identifying 
process P3 and determining an intersecting point and an 
intersecting line between them, it is possible to 
obtain geometric definitions of the vertex and edge 

25 line. 

Next, a construction of the X-ray CT 
apparatus used in the preliminary imaging process P2 is 
shown in Fig. 3. This X-ray CT apparatus has an X-ray 
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source (radiation source) Al, a scanning mechanism A2, 
a turntable A3, an X-ray detector A4, a control and 
signal processor A5, a marker arrangement pattern 
presenting means A6, and a marker arrangement pattern 
5 selection means A8 . 

X-rays emitted from the X-ray source Al 
penetrate an object (not shown) fixed on the turntable 
A3 and reach the X-ray detector A4 where they are 
converted into a detection signal XI. The detection 

10 signal XI is sent from the X-ray detector A4 to the 
control and signal processor A5 where it is subjected 
to "reconstructive processing" to be converted into a 
tomographic image set (three-dimensional bit-map data) . 

The turntable A3 is rotated and translated by 

15 the scanning mechanism A2 according to a control signal 
X2 from the control and signal processor A5- The 
rotation is, for example, a movement in a direction of 0 
in Fig. 3. The translation includes, for example, a 
movement in an X direction of Fig. 3 (lateral direction 

20 or longitudinal direction) and a movement in a Z 
direction (vertical direction) . By rotating or 
translating the turntable A3 with an object fixed on 
it, it is possible to radiate the X rays against the 
object from various directions and collect data 

25 (information) required for the reconstructive 
processing , 

Unlike the conventional X-ray CT apparatus, 
the X-ray CT apparatus of this embodiment has the 



- 13 - 

marker arrangement pattern presenting means A6 and the 
marker arrangement pattern selection means A8 . The 
marker arrangement pattern presenting means A6 has a 
marker arrangement pattern memory means A7 . The marker 
5 arrangement pattern memory means A7 stores a plurality 
of marker arrangement patterns in advance, such as 
shown in Fig. 2. The marker arrangement pattern 
presenting means A6 has a display means (not shown) 
such as a CRT display and displays a list of marker 

10 arrangement patterns on this display. The display 
means may also double as a display means used to 
control the X-ray CT apparatus. 

An example of screen the marker arrangement 
pattern presenting means A6 displays on the display 

15 means is shown in Fig. 4. On this screen are shown 
icons Vla-Vlc, a scrollbar V2, a selection button V3, 
and a cursor (arrow cursor) V4 . In Fig. 4, the icons 
Vla-Vlc correspond to the marker arrangement patterns 
shown in Fig. 2A-Fig. 2C- . 

20 The marker arrangement pattern selection 

means A8 comprises a general input means such as mouse 
and keyboard- The user may move the cursor V4 on the 
screen using the mouse to drag the scrollbar V2 to 
display a desired marker arrangement pattern on the 

25 screen. Then, the user selects an icon corresponding 
to the desired marker arrangement pattern with the 
cursor V4 and then clicks on the selection button V3 to 
specify the desired marker arrangement pattern. The 
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user can also specify a desired marker arrangement 
pattern by entering a number corresponding to the 
desired marker arrangement pattern from the keyboard. 

Because of the marker arrangement pattern 
5 presenting means A6, this embodiment can present the 
user with an appropriate marker attaching method. The 
provision of the marker arrangement pattern selection 
means A8 allows the user to select an appropriate 
marker attaching method easily. 

10 While in this embodiment we have shown an 

example case in which the CT apparatus has both the 
marker arrangement pattern presenting means A6 and the 
marker arrangement pattern selection means A8, the CT 
apparatus may have only one of them if it has a means 

15 for specifying a marker arrangement pattern. 

Next, details of the marker identifying 
process P3 will be explained by referring to Fig. 5. 
The marker identifying process P3 includes a density 
threshold setting process P31, a marker extraction 

20 process P32, a marker surface voxel selection process 
P33, and a surface domain expansion process P34. 

The density threshold setting process P31 is 
a process to set appropriate upper and lower thresholds 
for the marker density value. For example, when a 

25 plastic clay marker is to be attached to an aluminum 

object to be imaged^ a threshold for lower limit is set 
at 1.0 g/cm^ and an threshold for an upper limit is set 
at 2.0 g/cm^, taking into account the fact that the clay 



density is about 1-5 g/cm^. 

The marker extraction process P32 is a 
process which involves making comparison, for all the 
voxels making up the coarse tomographic image set Dl, 
5 between a material density described by the voxels and 
the threshold set by the density threshold setting 
process P31 and then extracting only those voxels 
corresponding to the markers. The voxel is a unit that 
makes up three-dimensional bit-map data. 

10 In the above example, when the densities 

described in the voxels are higher than the upper . 
threshold of 2.0 g/cm^, they are replaced with 0 g/cm^. 
This can eliminate the aluminum portion (voxels) from 
the coarse tomographic image set Dl, leaving only the 

15 marker portion (voxels) in the set- This process 
produces a marker tomographic image set D2 . 

The marker surface voxel selection process 
P33 is a process to pick up some voxels from those on 
the surface of the marker. First, this process 

20 randomly extracts voxels making up the marker 

tomographic image set D2 and checks whether the 
extracted voxels are on the marker surface or not. If 
the extracted voxels are on the surface of the marker, 
these voxels are added to a voxel list D7 . Here, 

25 several tens or several thousands of voxels are picked 
up - 

The surface domain expansion process P34 
assumes that a part of the marker surface has a 
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geometric feature, and expands a feature surface domain 
by taking one voxel included in that part as a starting 
point, extracts a domain satisfying the geometric 
feature, and determines the geometric feature. This 
5 processing is performed by using the surface domain 
expansion method. The surface domain expansion method 
can extract domains with geometric features (geometric 
definitions), such as a part of flat surface, a part of 
cylindrical surface, a part of spherical surface, a 

10 part of revolution surface, a part of quadratic 
surface, and a part of translated surface. 

Since the voxels on the marker surface are 
already extracted in the marker surface voxel selection 
process P33, these may be used as starting points . 

15 With the surface domain expansion method, it is 
necessary to determine beforehand which of the 
geometric features described above the domains 
containing the voxels of concern have, or to determine 
if these domains do not have such geometric features at 

20 all. 

For this purpose, the surface domain 
expansion process P34 first assumes that . a certain 
domain is a part of a planar surface and checks if this 
assumption is correct or not. This is done, for 
25 example, by attempting to extract a parameter defining 
the planar surface by the surface domain expansion 
method and then checking whether or not the parameter 
extraction is successfully completed. If it is found 
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that the domain is not a part of the planar surface, 
another assumption is made that the domain is a part of 
a cylindrical surface and then a similar check is made 
to see if this assumption is correct- In this way, by 
5 successively repeating the process of making an 
assumption and checking the validity of this 
assumption, it is possible to determine what geometric 
feature a domain including the voxels of concern has. 

Here, referring to Fig. 6A to Fig. 6D, one 

10 example process of extracting a geometric feature using 
the surface domain expansion method will be explained. 
In this example, a certain domain is assumed to be a 
"planar surface with an unknown normal." First, the 
coarse tomographic image set (three-dimensional bit-map 

15 data) Dl is transformed into point data Dla. An 
example transformation performed by a method that 
determines centers of two adjoining cells with a marker 
surface therebetween is shown in Figs. 6A-6D- Figs. 
6A-6D is a partially enlarged view of the coarse 

20 tomographic image set Dl, with cells (voxels) outside 
the marker shown blank and cells inside the marker 
shown hatched. Blank and black circles (O and •) at 
the cell boundaries correspond to the point data Dla. 

In this example, a group of point data Dla 

25 belonging to the "planar surface with an unknown 
normal" is determined and then a surface that best 
represents a set of points included in this group 
(representative surface) is determined. First, a group 
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addition process starts the following processing by- 
focusing on one cell. In the case of Fig. 6A, a cell 
bl is the one under consideration. A blank circle on 
the boundary of the cell bl is taken as a starting 
5 point SI, and two points adjoining SI are added to the 
group to expand the group. A state in which the group 
has been expanded to three points is shown in Fig. 6A. 
In Figs. 6A~6D, points belonging to the group are shown 
by "O" and those not belonging to the group by 

10 Next, a group boundary determination process 

determines a representative surface that best 
approximates a set of points belonging to the group by 
using the method of least squares. The representative 
surface is then offset toward the front and back 

15 (vertically up and down in Fig. 6B) to determine a 
boundary of the group. The amount of offset is 0.5- 
0-75 times a width of the cell. Fig. 6B shows a 
representative surface of the group indicated by a 
thick line N and its boundaries by thin lines Nl. 

20 Next, a group expansion process adds new 

points to the group based on the group boundaries Nl. 
More specifically, from among points not belonging to 
the group, those points adjoining the group are chosen 
and their coordinates are checked against the group 

25 boundaries Nl. If there is any points positioned 

inside the boundaries Nl, these points are added to the 
group. Fig. 6C shows an expanded state of the group 
after new points have been added to the group based on 



the boundaries Nl. 

After new points have been added by the group 
expansion process, the processing returns to the group 
boundary determination process where it recalculates 
5 the representative surface and the boundaries of the 
group that have changed because of the addition of new 
points. Fig. 6D shows a recalculated representative 
surface N and boundaries NX of the group. The group 
boundary determination process and the group expansion 

10 process are repeated until new points are no longer 
added to the group . 

When the points belonging to the group no 
longer change, a reference point and a normal of the 
representative surface obtained at that point in time 

15 constitute a geometric feature of this group. The 

reference point of the representative surface may, for 
example, be a point among a set of points belonging to 
the group that is situated at the center. This 
reference point and the normal are parameters, defining 

20 the planar surface. 

In determining the representative surface in 
the group boundary determination process, if a residual 
obtained as a result of the minimum square method fails 
to converge to a sufficiently small value, the 

25 assumption that the domain under consideration is a 
part of a planar surface is not valid. This is 
equivalent to a parameter extraction failure. 

Also for other point data Dla on the marker 



surface, the group and its representative surface are 
calculated in the similar manner to clarify the contact 
surface between the marker and the object. 

Described abbve is an example of how a 
5 geometric feature of the "planar surface with an 
unknown normal" is determined by the surface domain 
expansion method. The surface domain expansion method 
can determine a geometric feature according to a kind 
of the domain under consideration in a similar manner 

10 described above. 

For example, as a geometric feature of a 
"cylindrical surface with an unknown radius," a 
reference point, axis and radius are determined. For a 
"spherical surface with an unknown radius" a center and 

15 a radius are determined as its geometric feature. For 
a "quadratic surface" coefficients of a quadratic 
equation with X, Y and Z as unknowns are determined as 
its geometric feature. For a "revolution surface" a 
reference point of a rotating axis, a reference vector 

20 and a line being revolved are determined as its 
geometric feature. For a "translated surface" a 
translated vector and line are determined as its 
geometric feature. These geometric features also 
constitute information representing a geometry of the 

25 contact surface between the marker and the object. 

The surface domain expansion process P34 
produces a marker surface feature domain list D3, a 
list of domains on the marker surface having geometric 
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features. The marker surface feature domain list D3 
lists up domains of the marker surface that are larger 
than a predetermined area and which have certain 
geometric features. Of the domains on the marker 
5 surface, those not in contact with the object are 

generally free curved surfaces not classified into any 
of the planar surface, cylindrical surface, spherical 
surface, revolution surface and translated surface, and 
therefore are not recorded in the marker surface 

10 feature domain list D3 . 

A typical example of the marker surface 
feature domain list D3 is shown in Fig. 7, The marker 
surface feature domain list D3 is a set of marker 
surface feature domain records D31 describing 

15 individual domains having geometric features. The 
individual marker surface feature domain records D31 
have a marker entry D32, a geometric definition entry 
D33, an area entry D34, and a center of gravity entry 
D35- These entries describe a marker to which the 

20 domain of concern belongs, a geometric definition of 
the domain, an area of the domain and a gravity center 
of the domain, respectively. The content of the entry 
is calculated by the surface domain expansion process 
P34 and recorded in the corresponding entry. 

25 Next, details of the coordinate system 

determination process P4 will be explained by referring 
to Fig. 8. A three-dimensional Cartesian coordinate 
system can be defined by determining a reference point 
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and three orthogonal reference axes\ For this purpose 
six parameters need to be specified. These six 
parameters are expressed as (x, z, (j), \\f) . Of 
these, (x, y, z) are parameters for determining a 
5 position of the reference point, and G, \\f) are 
parameters for determining an affine transformation 
matrix that specifies directions of the three reference 
axes. {^f Qf \|/) are generally known as an Eulerian 
angle. 

10 The coordinate system determination process 

P4 determines six parameters from the marker surface 
feature domain list D3 according to the marker 
arrangement pattern the user has selected using the 
marker arrangement pattern selection means A8 . For 

15 example, in the case of the marker arrangement pattern 
shown in Fig. 2A, an intersecting point of the three 
planes recorded in the marker surface feature domain 
list D3 is taken as a reference, point and, based on the 
coordinates of the reference point, the parameters (x, 

20 y, z) are determined. The parameters ((|), G, V|/) are 
determined by using a multidimensional convergence 
calculation technique, called a valley descending 
technique or hill climbing technique. For example, 
they can be determined such that a value of an 

25 evaluation function shown in Equation 1 is minimum. 



[Equation 1] 
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f{(j>,e,il/) = Min 



{ZAX'f +(ZBry -\-(ZCZ')\ 
{ZAX'f-^iZBTf +{ZCTf, 
(ZAry -^{ZBX'f -h{ZCZ')\ 
{ZATf +{ZBZ'f+{ZCX'f, 
{ZATf +{ZBX'f+{ZCrf, 
{ZAZ'f +{ZBTf +{ZCX'f, 



(Equation 1) 



In Equation 1, A, B and C represent three 
planes recorded in the marker surface feature domain 
list D3- X', Y' and Z' represent three planes that are 
5 obtained through transforming the three reference 
planes (XY plane^ YZ plane, ZX plane) preset with 
respect to the turntable A3 of the X-ray CT apparatus 
by an affine transformation matrix determined by the 
parameters ((j), 9, xj/) . Z represents an angle between 

10 the intersecting planes . 

In this embodiment, a prespecified reference 
plane of an object is extracted by a CT image and then, 
based on a positional relation between this reference 
plane and the coordinate system of the turntable of the 

15 CT apparatus, a tomographic image of a prespecified 
portion is obtained. 

In the first embodiment described above, even 
if objects with different reference planes are set on 
the turntable for imaging, the user can omit the step 

20 of determining reference planes of each object, 

facilitating the placement of the objects and the 
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operation of the X-ray CT apparatus - 

In this embodiment an example case has been 
described in which the marker density is smaller than 
that of the object- Considering an attenuation of X- 
5 rays by the marker, this configuration is preferred. 
But if the density of the maker is larger than that of 
the object, this embodiment is still effective in 
omitting the work of determining the reference planes 
of the object. 

10 {Second Embodiment) 

In the first embodiment, we have described a 
case where plastic clay markers are directly attached 
to an object. A second embodiment concerns a case 
where markers are attached to an object through a pad. 

15 The pad has a small thickness compared with its size 
(length and width) and is made of a material with a 
smaller density than that of the marker. 

How the marker is attached to an object in 
this embodiment is shown in Fig. 9. The marker M 

20 includes a marker body Ml and a marker pad 1X12. The 

marker body Ml is made of a plastic clay with a density 
of about 1.5 g/cm^ and a size of 20-50 mm. The marker 
pad M2 is made of a light and deformable material such 
as foamed urethane, and is about 2 mm thick. 

25 The marker pad M2 has so small a density of 

about 0.1 g/cm^ that it does not show up in a tomogram 
output from the X-ray CT apparatus. Thus, on the 
tomographic image set, the object T and the marker body 
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Ml are completely separated, facilitating the. 
subsequent image processing to a great extent. To 
prevent the marker pad M2 from being imaged in the 
tomogram, the density of the marker pad M2 need only be 
5 20% or less of the density of the object T. 

The second embodiment of the CT imaging 
method using the marker body Ml and the marker pad M2 
of Fig. 9 is shown in Fig. 10. A point in which this 
CT imaging method differs from the first embodiment of 
.10 Fig. 1 is that this method includes a marker 
identifying process P3a instead of the marker 
identifying process P3. Other processes are the same 
as in Fig. 1 and are not described here. 

Details of the marker identifying process P3a 

15 will be explained by referring to Fig. 11. What is 
different from the marker identifying process P3 of 
Fig. 5 is that the marker identifying process P3a 
includes an opposing domain extraction process P34a 
instead of the marker surface voxel selection process 

20 P33 and the surface domain expansion process P34. 

Other processes are the same as those of Fig. 5 and are 
not described here. 

The opposing domain extraction process P34a 
first extracts domains of the surface of the marker 

25 body Ml whose distances from the surface of the object 
T are almost equal to the thickness of the marker pad 
M2. Next, domains of the surface of the object T that 
oppose the extracted surface domains of the marker body 



- 26 - 

Ml with the marker pad M2 interposed therebetween are 
determined. Based on the surface domains of the object 
T thus determined, the marker surface feature domain 
list D3, a list of domains having geometric features, 
5 is prepared. This processing can be performed, as in 
the surface domain expansion process P34 of Fig. 5, by 
successively checking whether or not the hypothesis 
that the domain of concern has a certain geometric 
feature is valid, 

10 As in the first embodiment, when objects with 

different reference planes are set on the turntable for 
imaging, this embodiment also allows the user to omit 
the step of determining the reference planes of the 
objects, facilitating the placement of the objects and 

15 the operation of the X-ray CT apparatus - 
(Third Embodiment) 

A third embodiment of the present invention 
concerning a CT imaging service method using the CT 
imaging method is shown in Fig. 12. This CT imaging 

20 service is for a service provider S to take tomographic 
images of an object T owned by a customer U by using an 
X-ray CT apparatus A and then provide the customer U 
with obtained tomographic images D5 and/or a 
tomographic image set D6, With this embodiment, the 

25 use of the method shown in Fig. 1 or Fig. 10 can not 
only simplify the tomographic imaging process but 
provide the CT imaging service also for objects that 
require a high level of confidentiality- 



- 27 - 

In this CT imaging service, a business 
transaction is initiated by a customer U applying to 
the service provider S for the service. Upon receipt 
of an application from the customer, the service 
5 provider S issues an instruction D8 to the customer U 
and sends a marker M or a plastic clay, the material of 
the marker, to the customer. The instruction D8 is a 
documentation that describes how to attach the maker M 
to an object T. The instruction D8 illustrates 
10 multiple marker arrangement patterns as shown in Fig. 
2. 

The customer U chooses an appropriate marker 
arrangement pattern from those illustrated in the 
instruction D8 and, according to the chosen pattern, 

15 attaches the marker M to the reference surfaces of the 
object T- That is, the customer U executes the marker 
application process PI shown in Fig. 1 and Fig. 10. 

After having attached all the necessary 
markers M to the object T, the customer U packs it in a 

20 shock-absorbing material, puts it in a container and 
sends it to the service provider S . At this time, the 
customer U prepares an instruction D9 specifying 
portions to be imaged and sends it along with the 
object T to the service provider S. The imaging 

25 instruction D9 includes information specifying the 

marker arrangement pattern selected by the customer U 
and information for specifying portions of the object 
the customer U wants to image. 
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The service provider S receives a container 
with the object T in it and fixes the container as is 
on the turntable A3 of the X-ray CT apparatus A. The 
container may be hermetically closed or sealed so that 
5 the object T may not be seen from outside. 

Next, the service provider according to 
the imaging instruction D9 from the customer U,^ enters 
information on the portions to be imaged into the X-ray 
CT apparatus A. The marker arrangement pattern 
10 specified in the instruction D9 is also entered into 
the X-ray CT apparatus A by the marker arrangement 
pattern selection means A8 . With these preparations 
completed, the image by the X-ray CT apparatus A is 
started- 

15 The X-ray CT apparatus automatically executes 

all the processes, i.e., the preliminary imaging 
process P2, marker identifying process P3 or P3a, 
coordinate system determination process P4 and real or 
final imaging process P5 . After the X-ray CT apparatus 

20 A has output the tomographic images D5 and/or 

tomographic image set D6 and their contents have been 
confirmed, the service provider S takes the container 
off the turntable A3 and returns it to the customer. 
The tomographic images D5 and/or tomographic image set 

25 D6 produced by the image operation are also returned to 
the customer U and a service fee is charged to the 
customer. The tomographic images D5 and/or tomographic 
image set D6 to be sent to the customer U are specified 



by the customer U at time of application. 

This embodiment extracts by the CT image the 
reference planes of an object specified by the customer 
U and, based on a relation between the reference planes 
5 and the coordinate system of the turntable of the CT 
apparatus, produces tomographic images of a portion 
specified by the customer 

In this CT imaging service, the object is 
returned in the same state as it was placed in the 

10 container by the customer. Thus, should a third party 
open the container and take the object during transport 
or storage, the customer U can easily detect an anomaly 
and take necessary actions appropriately. Further, if 
the container is sealed, the container cannot be opened 

15 by a third party without being readily noticed. This 
is suited to a CT imaging service for those objects 
that require high levels of confidentiality. Further, 
if the object has a complex shape, there is no need to 
manufacture a special jig for fixing the object at an 

20 appropriate angle, making it easier for the customer U 
to use the service . 

On the side of the service provider S, this 
service method saves the service provider the labor of 
picking up an object from the container and putting it 

25 again in the container. Further, this service method 
can eliminate a troublesome work of determining 
reference planes of the object required to align the 
angle of the object to the turntable. Therefore, the 



placement of the object and the operation of the X-ray 
CT apparatus are made easier, reducing the delivery 
period and service cost of the CT imaging service. 

In this embodiment an example case has been 
5 described in which a customer U prepares a. container to 
accommodate an object. It is also possible for the 
service provider S to send the container and shock- 
absorbing material to the customer U- In that case, 
the service provider S selects an appropriate size of 

10 the container and an appropriate material (density) of 
the shock-absorbing material according to information 
on the size and material (density) of an object 
obtained from the customer U. This permits the service 
provider S to efficiently execute the imaging of the 

15 object using the X-ray CT apparatus- 

The container suitably has a high X-ray 
transmissivity, a low X-ray attenuation rate, and a low 
density. More specifically, paper, corrugated 
cardboard, wood, plastics and resin may be used as 

20 container materials. A low density metal may also be 
used. 

In the above embodiments, we have described a 
CT apparatus using X-rays. It is also possible to 
apply the present invention to CT apparatuses using 
25 other radiations (a, P and y rays) and still produce the 
similar effects or advantages. In that case, the only 
requirement is to use a marker of a material with a 
different radiation transmissivity (or attenuation 
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rate) from that of an object to be imaged. 

With this invention, even when objects with 
different reference planes are to be set on the 
turntable for CT imaging, the user can omit the step of 
5 determining the reference planes of the objects and 
thus facilitate the placement of the objects. 

It should be further understood by those 
skilled in the art that although the foregoing . 
description has been made on embodiments of the 
10 invention, the invention is not limited thereto and 
various changes and modifications may be made without 
departing from the spirit of the invention and the 
scope of the appended claims. 



